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The interaction of a charged particle with its own field results in the "self- force" on the particle, 
which includes but is more general than the radiation reaction force. In the vicinity of the particle 
in curved spacetime, one may follow Dirac and split the retarded field of the particle into two parts, 
(1) the singular source field, ~ q/r, and (2) the regular remainder field. The singular source field 
exerts no force on the particle, and the self-force is entirely caused by the regular remainder. We 
describe an elementary multipole decomposition of the singular source field which allows for the 
calculation of the self-force on a scalar-charged particle orbiting a Schwarzschild black hole. 

O ■ I. INTRODUCTION 

(N 

Ph ' According to the equivalence principle in general relativity, a particle of infinitesimal mass orbits a black hole of 
^ large mass along a geodesic worldline T in the background spacetime determined by the large mass alone. For a 
particle of small but finite mass, the orbit is no longer a geodesic in the background of the large mass because the 
particle perturbs the spacetime geometry. This perturbation due to the presence of the smaller mass modifies the 
orbit of the particle from an original geodesic in the background. The difference of the actual orbit from a geodesic 
— 1| in the background is said to result from the interaction of the moving particle with its own gravitational field, which 
^ is called a self-force (author?) [Tl. 

y Historically, Dirac (author?) [2] first gave the analysis of the self- force for the electromagnetic field of a particle 

^ in flat spacetime. He was able to approach the problem in a perturbative scheme by allowing the particle's size to 
remain finite and invoking the conservation of the stress-energy tensor inside a narrow world tube surrounding the 
particle's worldline. Dewitt and Brehme (author?) [3] extended Dirac's problem to curved spacetime. Mino, Sasaki, 
i-H and Tanaka (author?) Q generalized it for the gravitational field self-force. Quinn and Wald (author?) and 
'i^ Quinn (author?) fo'l worked out similar schemes for the gravitational, electromagnetic, and scalar field self-forces by 
taking an axiomatic approach. 
' In Dirac's (author?) [2] flat spacetime problem, the retarded fleld is decomposed into two parts: (i) The flrst part 
is the "mean of the advanced and retarded fields" which is a solution of the inhomogeneous field equation resembling 
the Coulomb q/r piece of the scalar potential near the particle, {ii) The second part is a "radiation" field which is a 
homogeneous solution of Maxwell's equations. Dirac describes the self-force as the interaction of the particle with the 
radiation field, a well-defined vacuum field solution. 

In the analyses of the self-force in curved spacetime, the Hadamard form of Green's function (author?) [3] is 
employed to describe the retarded field of the particle. Traditionally, taking the scalar field case for example, the 
retarded Green's function G'^'^^{p,p') is divided into "direct" and "tail" parts: (i) The first part has support only on 
the past null cone of the field point p, {ii) The second part has the support inside the past null cone due to the 
^ ' presence of the curvature of spacetime. Accordingly, the self-force on the particle would consist of two pieces: (i) The 
5^ , first piece comes from the direct part of the field and the acceleration of the worldline in the background geometry; 
this corresponds to Abraham-Lorentz-Dirac (ALD) force in fiat spacetime, {ii) The second piece comes from the tail 
part of the field and is present in curved spacetime. Thus, the description of the self-force in curved spacetime should 
reduce to Dirac's result in the flat spacetime limit. In this approach, the self-force is considered to result via 
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= qVa^J, (1) 

from the interaction of the particle with the quantity (author?) [l| 

Although this traditional approach provides adequate methods to compute the self-force, it does not share the 
physical simplicity of Dirac's analysis where the force is described entirely in terms of an identifiable, vacuum solution 
of the field equations: unlike Dirac's radiation field, the quantity in Eq. (??) is not a homogeneous solution of the 
field equation V'^tp — —Airg. Moreover, the integral term in this quantity comes from the tail part of the Green's 
function and is generally not differentiable on the worldline if the Ricci scalar of the background is not zero (similarly, 
the electromagnetic potential ^Jj**'' and the gravitational metric perturbation /ijj^'' are not differentiable at the point 
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of the particle unless (i?a6 — ^9abR) and Rcadbu'^u'^ , respectively are zero in the background (author?) Q). Thus, 
some version of averaging process must be invoked to make sense of the self-force. 

In this paper an alternative method to split the retarded field ■0'^°* as suggested by Ref. (author?) fl\ is used 
such that the splits resemble those by Dirac: (i) The first part, the Singular Source field is an inhomogeneous 
field similar to the Coulomb potential piece and exerts no force on the particle, (ii) The second part, the Regular 
Remainder field ijj^, being an homogeneous solution of the field equation, is analogous to Dirac's radiation field and 
entirely responsible for the self- force. This alternative method is reviewed briefly in Section HIl 

In Section IIIII we give a brief overview of the mode-sum decomposition scheme to evaluate the self-force. We 
consider a particle with a scalar charge q in general geodesic motion about a Schwarzschild black hole in this paper. 
A spherical harmonic decomposition of both ip'^'^^ and tp^ having this condition should be performed to provide the 
multipole components of each. Then, the mode by mode sum of the difference of these components determines 
and, thence the self-force. The multipole components of ip'^'^^ can be determined numerically while the multipole 
components of are derived analytically. In particular, the multipole moments of tp^ are generically referred to as 
the Regularization Parameters; the entire paper focuses on the analytical task to find these regularization parameters. 
We summarize our analytical results at the end of the Section. 

The description of ilr" becomes advantageously simple in a specially chosen co-moving frame: the THZ normal 
coordinates of this frame, introduced by Thorne and Hartle (author?) Q and extended by Zhang (author?) Q 
are locally inertial on a geodesic. In Section ITVl we obtain a simple form of ip^ using the THZ coordinates and then 
re-express it in terms of the background, i.e. Schwarzschild coordinates via the coordinate transformation between 
the THZ and the background coordinates. 

Section |V] outlines the derivation of the regularization parameters given below in Eqs. (ITU) to (|18l) . These results 
are in agreement with Barack and Ori (author?) [10] and Mino, Nakano, and Sasaki (author?) [11]. 

In Appendix we give a detailed description of the THZ coordinates. Appendix [B] provides some mathematical 
details concerning the hypergeometric functions and the different representations of the regularization parameters in 
connection with them. 

Notation: (t, r, 9, (p) are the usual Schwarzschild coordinates. (T, X, Y, Z) are the initial static normal co- 
ordinates, reshaped out of the Schwarzschild coordinates. (T, X, y, Z) are the THZ normal coordinates attached 
to the geodesic F, and p = \/ X"^ +y'^ + Z'^. The points p and p' refer to a field point and a source point on the 
worldline of the particle, respectively. In the coincidence limit p ^ p' . 



II. DECOMPOSITION OF THE RETARDED FIELD 

The recent analysis of the Green's function decomposition by Detweiler and Whiting (author?) [1] shows an 
alternative way to split the retarded field into two parts 



7/,^<=t = V^s + (3) 

where ip^ and ^p^ are named the Singular Source field and the Regular Remainder field, respectively. The source func- 
tion for a point particle on the worldline T is g{p) — q J{—g)^^^'^5^{p — p'{T'))dT' . Like ^p^ is an inhomogeneous 
solution of the scalar field equation 



V^V = -Ang (4) 

in the neighborhood of the particle. And ip^ is determined in the neighborhood of the particle's worldline entirely 
by local analysis, ip^, defined by Eq. ([3]), is then necessarily a homogeneous solution and is therefore expected to be 
differentiable on F. According to Ref. (author?) ip^ will formally give the correct self-force when substituted on 
the right hand side of Eq. ([1]) in place of ■0*'^''. In this paper we adopt this decomposition, and try to determine an 
analytical approximation to ip^, which is to be subtracted from tp'-'°^ for an explicit computation of the self- force. 



III. MODE-SUM DECOMPOSITION AND REGULARIZATION PARAMETERS 



By Eq. ([T]) the self-force can be formally evaluated from 
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q lim VaV'^ = 9 lim Va(V''''* - V'^) 

p—^p' p—^p' 



(5) 



where p' is the event on F where the self- force is to be determined and p is an event in the neighborhood of p' . For 
use of this equation, both F'^^'^{p) and T^{p) would be expanded into multipole ^- modes, with J^l'^{p) determined 
numerically. 

Typically, if the background geometry is Schwarzschild spacetime, the source function q{p) is expanded in terms of 
spherical harmonics, and then similar expansion for if/''^ is made 



r 



(6) 



where ip^'^ (r, t) is found numerically. The individual ^m components of in this expansion are finite at the location 
of the particle even though their sum is singular. Then, J-^^ is finite and can be obtained as 



(7) 



where a represents each component of t, r, 0, 9 in the Schwarzschild geometry. 

The singular source field ij}^ is determined analytically in the neighborhood of the particle's worldline via local 
analysis (see Section ITV)) . Then, VaV'^ is evaluated and the mode-sum decomposition of this quantity is performed 
to provide 



im 1 



(8) 



which is also finite at the location of the particle. 

Then, using Eqs. ([5]), jT]), and ([8]) the self- force is finally 



P. 



■self 



5]lim^ [Flf{p)-Fl{p)\ 



(9) 



evaluated at the location of the particle. 

In Section |V] the regularization parameters are derived from the multipole components of Vqt/i^ evaluated at the 
source point and are used to control both singular behavior and differentiability. We follow Barack and Ori (author?) 
[l^ in defining the regularization counter terms, except that the singular source field i/;^ is used in place of ■0'^"^ 



and show 
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Aa+Ba + ^ + 0{r^), 



(10) 



At = sgn(A) 



r2 1 + Pjrl' 



(11) 



Ar = -sgn(A) 



2 ^ ( 1 - 2M 



rl 1 + Pjrl 



(12) 



= 0, 



(13) 
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5/2 



{l + J^/rlf^ 2(1 + j2/r2)5/2 



(14) 
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[1-2(1- ^) r^]F,;2 3 (1 - MI) PF,/, ^ 



'^o I (l + j2/^2)i/2 2(1 + J2/r2)'/' 2(1 + J2/r2)5/2 

-P'l/2 - 3(-p5/2 " ^3/2) 



(15) 



(16) 



(l + JVr2)i/2 2(1 + j2/r2)3/2 

Ct=Cr = C^ = 0, (17) 

= Be = = 0, (18) 

where A = r — Tq, = — = (1 — IMjro) [dt/dr)^ (r: proper time) and J = Ufj, = {dcji/dr)^ are the conserved 
energy and angular momentum, respectively, and f = = (dr/dr)^. Also, shorthand notations are used for the 

hypergeometric functions, Fp = 2F1 {p, ^; 1; ^2'^ 72 ) (see Appendix |B] for more details about the hypergeometric 

functions and the representations of the regularization parameters in connection with them). 

IV. DETERMINATION OF VIA THE THZ NORMAL COORDINATES 

It was mentioned earlier in Sections lU and Hill that -0^ is determined in the neighborhood of the particle's worldline 
entirely by local analysis. When analyzed by some special local coordinate system in which the background geometry 
looks as flat as possible, the scalar wave equation might take a simple form and ij:^ might look like a simple Coulomb 
potential piece. Detweiler, Messaritaki, and Whiting (author?) ^] (cited henceforth as Paper I) show 

^S = 5/p + 0(p2/7^3)^ (19) 



where p = V ^"2 + 3^2 _|_ 2^2 with A", y , Z being spatial components in that local inertial coordinate system and TZ rep- 
resents a length scale of the background geometry (the smallest of the radius of curvature, the scale of inhomogeneities, 
and time scale for changes in curvature along V). 

Describing this special coordinate system more precisely, first, it must be a normal coordinate system where on 
r, the metric and its first derivatives match the Minkowski metric, and the coordinate T measures the proper time. 
Normal coordinates for a geodesic, however, are not unique, and we use particular ones that were introduced by 
Thorne and Hartle (author?) [8,] and extended by Zhang (author?) [§| to describe the external multipole moments 
of a vacuum solution of the Einstein equations: namely, the THZ NORMAL COORDINATES. 

However, in order to derive the regularization parameters from the multipole components of VaV'^j P in Eq- (O 
must be expressed in terms of the coordinates of the original background, which is the Schwarzschild geometry in our 
problem. Then, this requires us to find the expressions of A", 3^, Z (THZ) in terms oit, r, 0, 6 (Schwarzschild): the 
task here is simply to find the coordinate transformation between two different geometries. 

Based on the idea from Weinberg (author?) [l^, one can achieve this coordinate transformation to the level of 
accuracy we desire for this particular problem of mode-sum regularization, by taking the following two steps basically: 

(i) Find an inertial Cartesian coordinates X"^ to redirect the Schwarzschild coordinates x"' by the Taylor's expansion 
around the location of the particle, x"; 

= X^ + M^aix^ - + Im^, nX ix' - x'oKx' - O + 0[{x - Xo)% (20) 

where we may choose X^ = Q and M^a = diag [M'^t, M^r, Af^^, M^e] for convenience (this choice will 
redirect and rescale the Schwarzschild coordinates as T = M'^t{t — to), X — AI^ r{r ~ ro), Y — M^^ ^{(p — (po), 
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(ii) Boost X"^ with u"^, the particle's four- velocity at p' as measured in this Cartesian frame, to obtain the final 
coordinates : 



+ 0[{x-Xo)\ 



(21) 



where 



7~| V V 7 

u —u —u —u 

SYM l + {u^ -\){u^ f/u'^ {u'^ ^\)u^u^ /u^ 

l + (u^-l)(u^)Vu' 



(22) 



with u2 = („X)2 + (yi')2 ^^zy (author?) [T|. 
According to Ref. (author?) [13], one can show out of Eq. (|2ip 



.,A'B' 



ah 



dx'^ dx^ 
^'^'+O[(x-xo)2], x'^^xl, 



so that 



dg 



A'B' 



dX^ 



- =0[(S-Xo)], 



(23) 



(24) 



with the choice of M^t = " ^) ^'^^ " ~ ^) ^^^"^ foSin6lo, M^e = -ro. Eqs. ^ and ([241) 
are the local inertial features as expected for normal coordinates. 

To simplify the calculations, we may confine the particle's orbits to the equatorial plane 6^ — 7r/2. Then, we have 



-1/2 



M^a = diag 



1 ' O 1 ' O 



(25) 



where / = ^1 — . Also, this constraint of the equatorial plane makes = 0. We can rewrite in terms of the 
Schwarzschild coordinates and the constants of motion. 



u = [u , u , u , u 



(26) 



where E = —ut = f (dt/dr)^ (r: proper time) and J = = {d(j)/dT)^ are the conserved energy and angular 
momentum, respectively, and f = = (dr/dr)^. From this it follows that — f~^E^ — 1 and we have 



A^'a 



f-y^E -r^'^r -J/r, 

l + P/U^^E + f) Jf/[r,{E + f^/^)] 
SYM l + J^/[rl{f-^/^E+l)] 

1 



(27) 



Now we are finally able to express in terms of the Schwarzschild coordinates. Using Eq. (pi]) one may write 

p2 = x'Xi = SijA'c^'dM^^.M^^ [{x- - xDix" - xt) + Tl,l [x^ - x^ix"- - x^ix" ~ xt)] 

+0[{x~x,)% (28) 

where /, J = 1, 2, 3. Then, using Eqs. and (P7)) . Eq. (pS)) can be eventually expressed as 



{E^-f){t~t,f 



2Er 



{t - to)ir - To) - 2EJ{t - t^){cj) - cj)^) 
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+ J (l + 7) - + - -o)(0 - ^o) + {rl + .m - 4>of +rl{e-l)" 

MEr, Mf 2E^ P\ ^ x MJr, , 

MEr , ,y 2(ro-M)EJ, 

J ri fri 



-ro Jr(0 - (j^of - roJH^ - 0o) - f ) ' + - Xo)^]. (29) 

Eq. (1^^ is substituted into Eq. ([T^ to determine -0^ in terms of the Schwarzschild coordinates and will serve 
significantly to derive the regularization parameters in the next section. 

In the above analysis we have ignored the term 0[{x — Xq)^] in Eq. (|2ip and its contribution to p^, which is 
0[{x — Xo)*] in Eqs. and (P^ . To the level of accuracy we desire for the mode-sum regularization in this 
paper, that is to say, to the determination of C-terms, it is not necessary to specify the term 0[{x — Xo)'^] in 
(hence not necessary to specify 0[{x — Xo)^] in the spatial THZ coordinates X, y, Z). Even without specifying 
the term 0[{x — Xo)^] in p^, one can prove that Co-terms in Eq. ([TU)l always vanish (see Subsection IV C[) . In fact, 
the clarification of 0[{x — Xo)^] for the THZ coordinates in Eq. (PT|) requires more involved analyses of coordinate 
transformations, which would be beyond the scope of this paper. Readers may refer to Appendix |X] for more detailed 
description of the THZ coordinates, specified up to the quartic order. 



V. REGULARIZATION PARAMETERS FOR A GENERAL ORBIT OF THE SCHWARZSCHILD 

GEOMETRY 

In Section [IVI we have seen that an approximation to -0^ is 

^A' = (7/p + 0(pV7^')• (30) 

Following Paper I (author?) Q, the regularization parameters can be determined from evaluating the multipole 
components of da{q/ p) {a = t, r, 9, (p for the Schwarzschild background) at the location of the source. The error, 
0{p^ /Ti?) in the above approximation is disregarded since it gives no contribution to VaV'^ ^-s we take the "coincidence 
limit", X ^ Xq , where x denotes a point in the vicinity of the particle and Xo the location of the particle in the 
Schwarzschild geometry. 

In evaluating the multipole components of da{q/p), singularities are expected with certain terms. To help identify 
those singularities, we introduce an order parameter e which is to be set to unity at the end of a calculation: we attach 
e" to each 0[{x — Xo)"] part of p^ in Eq. and may re-express p^ as 

p2 - e^Vu + e^T'iii + e^Vw + O(e'), (31) 

where Vu, Vm, and Vrj represent the quadratic, cubic, and quartic order parts of p^, respectively. Here we pretend 
that the quartic part Viy is also specified: this will help us to perform the structure analysis for Ca-terms later in 
Subsection IV CI when we prove that these regularization parameters always vanish. 

When we express da i^/ p) in the Laurent series expansion to identify the terms according to their singular patterns, 

every denominator of this expansion should take the form of Vl\^^ (n = 3, 5, 7). Due to this special position, which 
becomes singular in the coincidence limit, Vu would play a significant role in inducing the multipole decomposition. 

However, the quadratic part Vu, directly taken from Eq. ((29|) . may not be fully ready for this task yet. First, (t) — (t>o 
needs to be decoupled from r — so that we have independent complete square forms of each, which is necessary for 
inducing the Legendre polynomial expansions later. Coupling between t — to and (j) — (pa is not significant because 
upon fixing t = all terms having t — to will vanish. Thus, we reshape our quadratic term in Eq. (|29l) into 

2Err^ 

Vu = {E^^m-tof- ° {t-to)^-2EJ{t~to){(l>-<l>') 
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+ - 



P {rl + J2) 
with 



^^ + {rl + j')[c^-<l>'f+rl{e^l)" (32) 



(33) 



/(r2 + j2)' 



where A = r — Tq, and an identity P = — f {l + J"^ is used for simphfying the coefficient of A^. Here, taking 
the coincidence hmit A 0, we have 0' — >■ 0o (the same idea is found in Mino, Nakano, and Sasaki (author?) [Hj). 
Also, in order to be muhipole-decomposed globally, the quadratic part needs to be completely analytic and smooth 
over the entire two-sphere. For this purpose we rewrite it as 



Pii 



2Efrl 



{t - to) A - 2EJ(t - to) sin sin(0 - (p') 



P {rl + J2) 
-0[{x-xof] 



A2 + (r2 + J2)sin' 



2^sin2( 



b') + rl cos^ 



(34) 



sm 



0') + O[(0 - cp'f] and 1 = sin 9 + 0[{e - 7r/2)2] we create 



where one should notice that by replacing cf) — < 
the 0\{x — XoY] terms. 

To aid in the multipole decomposition we rotate the usual Schwarzschild coordinates by following the approach of 
Barack and Ori (author?) [l^l and Paper I (author?) Q such that the coordinate location of the particle is moved 
from the equatorial plane = ^ io a location where sin0 = (9: new polar angle). We define new angles 9 and $ 
in terms of the usual Schwarzschild angles by 



sin6'cos(0 — (j)') 
sin0sin(0 — 0') 



cos 9 — 

Also, under this coordinate rotation, a spherical harmonic 



cos 9 

sin 9 cos $ 
sin 9 sin $. 



(35) 



becomes 



E 



,,y,„K9,<j>), 



(36) 



where the coefficients Q!f,„„/ depend on the rotation 



(9, $) as well as on ^, m, and to', and the index t 



is preserved under the rotation (author?) [l^j. As recognized already in Paper I (author?) [2,], there is a great 
advantage of using the rotated angles (9, $): after expanding da{q/ p) into a sum of spherical harmonic components, 
we take the coincidence limit A — > 0, 9 — 0. Then, finally only the m = components contribute to the self-force at 
9 = since Y«m(0, $) = for m 7^ 0. Thus, the regularization parameters of Eq. (fTO)) are just (^, to = 0) spherical 
harmonic components of da{q/p) evaluated at x^- 

Now, using these rotated angles, we may re-express Eq. (p4l) as 



{E^-f){t-tof 



2ErP 

^ (t - io) A - 2EJ{t ~ to) sin 9 cos $ 



+2 (P + P) 



+0[{x - Xo)' 



1 - 



fiP + P) 
P svP $ 



2/2 (r2 + J2)2 (1 _ J^^) 



1 — cos 9 



(37) 



where an approximation sin2 9 = 2(1 — cos 9) -I- 0(9**) is used, the error from which is essentially 0(9*^) = 0[{x — Xo)'^] 
and can be absorbed into Viy. Here one should note that through a series of modifications of the quadratic part of 
Eq. (|29p we have created additional quartic order terms apart from the desired form. Then, we may remove these 
additional terms from the quadratic part and incorporate them into the quartic part Viy in Eq. (|3ip . It is not 
necessary, however, to specify this quartic part: as already mentioned above, later in Subsection IV CI we will show 
that the quartic part Pjv starts appearing from the e°-term and verify that the regularization parameters of the 
e'^-term always vanish by analyzing the generic structure of the e'^-term. 
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After removing 0[{x — Xo)"*] from Eq. p7p we may define 

P' 

+2 (rl + j') I 1 



"2 - (E^ - f)(t - tof ~ —— — ^{t~to)A-2EJ{t-to)smQcos<^> 



J2 sin^ $ 



-2 + J2 



2/2 (r2 + ./2)2 _ i^^j 



1 — COS 6 



In particular, when fixing t — to, we define 



~2 _ -2 1 



2(7-0 + + i -cose) 



(38) 



(39) 



with 



J2 sin2 $ , , 

X = 1 - 2^ ,2 40 

and 

r2 R2 A2 
2/2(r2 + j2)2^ 

Now we rewrite Eq. (|3ip by replacing the original quadratic part Pn with the modified form above, 

^e^p^+e^Vui + e^Viv + 0{e''), (42) 

where the new quartic part Vw includes the additional quartic order terms that result from modification of the 
quadratic part Vu- Then, based on this redefined p^, we have the following expression of da{^/ p)\t=t ^ Laurent 
series expansion 



da ( - 
P 



2 1 2 p3 +4 ^5 r 



Oie'). (43) 



Eq. (|39p . when inserted into Eq. (|43p . plays a very important role in calculating the regularization parameters for 
the rest of the section: out of Eq. p9p . we induce Legendre polynomial expansions in terms of cosO. Sometimes we 
may have the dependence on not only in the denominators but also in the numerators on the right hand side of 
Eq. (j43| . In the numerators the dependence can be found from the terms containing sin" 8 or cos" Q since Eq. ((39|) 
can be solved for cosO. After finding all of the p^ dependence, the rest of the calculations involve integrating over the 
angle $ . The techniques involved in Legendre polynomial expansions and integration over <I> are described in detail 
in Appendices C and D of Paper I (author?) [7|. 

Below in Subsections IV Al and IV B[ we present the key steps of calculating the regularization parameters. 

A. Aa-terms 

We take the e^2 term from Eq. ((43|) and define 

Then, we proceed with our calculations of the regularization parameters one component at a time. 

1. At-term: 

First we complete the expression for Qt[e~2] by recalling Eqs. p8|) and (|39p 
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(7^ r / 9 9\ \ 1-3/2 / 2£'fr^A 

1- [2 {rl + J') xiS' + l -cose)] ^ ' 



f{rl + J^) 



2i?Jsin9cos$ 



((52 + 1 -cose)"'/'. 



A 



where means that A is held constant while the differntiation is performed with respect to Q. 

According to Appendix D of Paper I (author?) f?!, for p > 1 



(^2 _^ 1 _ cosQ)-"-"^ = J2 s2p^i(V I) [1 + PdcosQ), 6 ^ 



21+1 



and for p = Q 



{6^ + l-cosey^^^ [V2 + 0(M) P<;(cose), 5-^0. 

i=0 

Then, by Eqs. ge]) for p = 1, (gT]), and ([41]), in the Hmit 5^0 (equivalently A ^ 0) Eq. (gS]) becomes 



(45) 



(46) 



(47) 



lim Qt[e^^] = sgn(A) 

A->0 



2 • —1 / -, 



'■EJx-^/^ COS'S? ^ d 



(^2 + j2)3/2 99 



i=0 



PeicosQ). 



A 



(48) 



Then, we integrate limA->o Qt[e ^1 over $ and divide it by 2tt (we denote this process by the angle brackets "()") 



lim^Q,[e-2]) = sgniA)^p^^^ Y.(e+\] F,(cose), 



(49) 



where we exploit the fact that (x^^^^ cos $) = to get rid of the second part in Eq. (^51) fis']. Appendix C of Paper I 

(author?) [7] provides (x~') = 2F1 (l, ^; 1; a) = Fi — {I — a)"^^^, where a = J^/ (r^ + J^). Substituting this into 
Eq. (j49p , the regularization parameter At can be finally determined when we take the coefficient of the sum on the 
right hand side in the coincidence limit B 



At = sgn(A) 



r2 1 + J^/r"^ 



(50) 



2. Ar-term: 



Similarly, we have 



Here, before computing dr we reverse the process via Eqs. ([5^ . ([M)) . (|37p . and 



(51) 

to obtain the relation 



Vii + 0[{x-x,Y] 
E^rl 



2Efri 



(t~to)A-2EJ{t~to){cl)- 
2 



+O[{x-Xo)^] 



A' + {rl + J2) (0 - +rl{e- I) 



(52) 



(53) 
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Differentiating this with respect to r and going through the process via Eqs. p4p and psp . Eq. (|51l) can be expressed 
with the help of Eq. as 



^[2(r^ + j2)x(<52 + l-cose) 



-3/2 



J2 



/Jrsin0cos$ 



(54) 



[igj . Then, the rest of the calculation is carried out in the same fashion as for the case of Af-term above. We obtain 

„2 E 



-sgn(A) 



rini+j^/rir 



(55) 



3. AA-term: 



First we have 



Taking the same steps as used for A^-term above via Eqs. ([55)1 . and in order, we obtain 

^0 = 2 (r^ + J2) sinecos$ + 0[{x - x,f]. 

Then, in a similar manner to that employed in the previous cases, in the limit A ^ Eq. (|56p becomes 



lim Q^[e 



3/2 (,os$ ^ d 



PiicosQ) 



[20j . The right hand side vanishes through "()" process because cos$) = 0. Hence, 

= 0. 



(56) 



(57) 



(58) 



(59) 



4- Ag-term: 

It is evident from the particle's motion, which is confined to the equatorial plane 0o = f , that no self force is 
acting on the particle in the direction perpendicular to this plane. This is due to the fact that both the derivatives of 
retarded field and the singular source field with respect to 9 tend to zero in the coincidence limit. Our calculation of 
Ag should support this. Through the same process as employed before, we have 



with 



de (p^) 1 = sin 9 sin $ + 0[{x - x^f 



Then, similarly as in the case of Aw,-term above 



lim Qi}\t 

A->0 



(j^r\x ^^^sin$^ d 



Pf(cose). 



Again, via "()" process, the right hand side vanishes because (x^"^^^ sin$) = 0. Thus, 

A« = 0. 



(60) 



(61) 



(62) 



(63) 
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B. Ba-terms 

We take the e^^ term from Eq. ((43|) and define 



Qa[e ]=ql--^^ + - , (64) 



where for computmg da (p^) , Eq. ([53| should be referred to, and Vui is the cubic part taken directly from Eq. ([29]) . 
We may express this in a generic form 

2 2n [fc/2] , ^.2n-k/l A \'^-'^P ( t^\^P 

Q^i^''^ = E E E — — (65) 

n=l k=0 p=0 P° 

where A = r — Tq, and bnkp(a) is the coefficient of each individual term that depends on n, k and p as well as a, with 
a dimension 7?.''"^^ for a = t, r and 7?.''' for a — 6, (f). We recall from Eqs. (|32p and ([33| that the first of the steps to 
lead to in the denominator is replacing (j) — (f>o hy {(f> — <j)') ~ j(jd^rpj^ to eliminate the coupling term A (0 — 0o)- 
This makes a sum of independent square forms of each of A and (j) — 4>', which is a necessary step to induce the 
Legendre polynomial expansions later. Thus, to be consistent with this modification in the denominator, (j) — (f>o in the 
numerator on the right hand side of Eq. (|65|) should be also replaced by {(j) — (j)') — /(r^'+j^) Then, this will create 

a number of additional terms apart from (</> — (j)')™' when we expand the quantity {(f> — <j)') — ^^^^^p^ /S. raised, say, 

to the m-th power, and the computation will be very complicated. 

By analyzing the structure of the quantity on the right hand side of Eq. (|65p one can prove that — 0o may be 
replaced just by — 0' in the numerator without the term — j^^^^-j^jA (the same idea is found in Mino, Nakano, and 
Sasaki (author?) [11]). The verification follows. The behavior of the quantity on the right hand side of Eq. (I55t . 
according to the powers of each factor, is 

Qa[e-'] - p-t^^+DA^-'^ (0 - 0o)'-'^ (e - |)''7^^ (66) 
where s — k — 1 for a — t, r and s — k for a = 6, (j). Further, 

fe-2p 
1=0 

- (sine)* (cos $)''A'=-2p-Y7^fc-2p-» + o[(.T-Xo)'=-'''+'], (68) 

where a binomial expansion over the index i — 0, 1, • • • , fc — 2p is assumed with Ckpi ^ I/TZ'^^'^p^^ in Eq. (|67p . and in 
Eq. dSSl) ((/> - 0')' is replaced by [sm{(j) - </>')]' + O[(0 - 4)')'^'^] — the term 0[{x - Xo)''^'^^^'^] at the end results from 
this 0[{(p — (f>'y^'^], then the coordinates are rotated using the definition of new angles by Eq. (|55]) . Also, by Eq. ([55]) 
again 



= {smQfP (sm^)^" + 0[{x - Xo)^P+^]. (69) 

Using Eqs. ^ and ((Ml), the behavior of (3[e"^] in Eq. ([66]) looks hke 

Qa[e-^] - p-(2"+i)A2»~2p-^ (sine)'P+* (cos$)* (sin$)'P7^^ (70) 



where s = 2p + i — l for a = t, r and s = 2p + i for a = 6, (f), and any contributions from 0[(x — Xo)*'~^^^^] in Eq. 
and from 0[{x — Xo)^^'''^] in Eq. (IM]) have been disregarded: by putting these pieces into Eq. (I66p we simply obtain 
e^-terms, which would correspond to 0{i^^) in Eq. (llOp and should vanish when summed over £ in our final self- force 
calculation by Eq. ^ [2l|. Qa[e~^] then can be categorized into the following cases: 
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(i) i = 2j + l (j = 0, 1, 2, • • •) 

The integrand for "()" process, = (cos$)^^+^ (sin$)^'' has the property + n) = -F{^). Thus 

{Q4e-'])=0, (71) 

(ii) i = 2j {j = 0, 1, 2, • • •) 

Using Eqs. (|39|) and (j4T|) . we can express (sinG)^^"''' in Eq. (|70|) above in terms of po and A via a binomial 
expansion 

(sin 6)'^+'^' = [2 (1 - cos 6)]^+^' + 0[{x - x,f^P+^^+^] 

= E d„,pl^A''^P+^-^^ + 0[ix - Xo)^(''+^)+2] (72) 

9=0 

where g = 0, 1, • • • , p + j is the index for a binomial expansion and dpjq ^ When Eq. (|73p is 

substituted into Eq. (|70|) . the contribution from 0[(x — Xo)^'^"''-'-'^^] can be disregarded since it would correspond 
to O(e^) again. Then, we have 

Qa[e"'] - (sin $)'P (cos $)'^' p-2(„-g)-1^2(„-g)-^s^ (74) 

where s = — 1 for a ^ t, r and s = for a ~ 9, (/), and we can guarantee that n — g > always since 
0<(Z<p + i= p+ |i,0<i<fc — 2p, and p < k < 2n. Then, Eq. ([74| can be subcategorized into the following 
two cases; 

(ii-1) n-q>l 

By Eqs. HMD, gU, and (gSl) 

Qa[e^^] - (sin$)^P(cos$)^^' AP^(cose)7^" ^ 0, (75) 

A-i-O 

(ii-2) n-q = 

By Eqs. HMD, 1111), and dUD 

Qa[e-i] - (sin$)'P(cos$)'^>^(cose)7^^ (76) 

where s = — 1 for a — t, r and s = for a ^ 9, (f). 

Therefore, by analyzing the structure of Qa[e^^] we find that the e~^-terms vanish in all the cases except when 
n — q = {). The non-vanishing _Ba-terms are derived only from this case. Then, by < q < p + j — p + ^i, 
< i < k — 2p, and p < k < 2n together with n = q one can show that 

< fc - 2p - i and fc - 2p - i < 0, i.e. k - 2p - i = 0. (77) 

Substituting this result into Eq. (p7)) . then into Eq. (|65l) we may conclude that in the numerator of Q[e^^] in Eq. ([55)1 
one can simply substitute 

(0 - 0o)'"'^ (0 - q^'f-'" ■ Q. E. D. (78) 

The significance of this proof does not lie in the result given by Eq. ((78| only, but also in the fact that the 
non- vanishing contribution comes only from the case n = q for Eq. (17411 . i.e. 

Qa[e-'] ~ ism<S>fP (cos $)"("-") p-'U' , (79) 

where n = 1, 2 and < p < n, and s = — 1 for a = t, r and s ~ for a ~ 9, (f). 

Below are presented the calculations of _Ba-terms of the regularization parameters by component, in a similar 
manner to those for A„ -terms . 
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1. Bt-term: 



We begin with 



(80) 



, 2 p3 4 

The subsequent computation wih be very lengthy and it will be reasonable to split Qf[£~^] i^ito two parts. First, let 

Qtii){e-']^-^Po' dtVui\t=t^. (81) 



where 



2 1 



M\ EJA 



fro 



TnE'r 



(82) 



As proved at the beginning of this Subsection, every (</> — 0o)™ in the numerators of the e ^-term can be replaced by 
((/) — f/)')™ without affecting the rest of calculation. Then, followed by the rotation of the coordinates via Eq. ([551 



Qt(i)[e" 







MErA 


2 




[ Prl 


+0 


Ux- 


-Xof] 


1 



2(1- — 



M\ EJA 



ho 



sin 9 cos $ + 2roEr (1 — cos 6) 



(83) 



where an approximation sin^ Q = 2(1 — cos Q) + 0[{x — Xo)'^] is used to obtain the last term inside the first bracket. 
Here we may drop off the term O [{x — Xo)^/p^'j , which is essentially O(e^), for the same reason as explained at the 
beginning of this subsection. Then, using the same techniques as used to find Aa-terms, we can reduce Eq. (|83p to 



Qt 



q^MEr 



q r^E rx 



V^rl 2/2 {rl + .py 



1 - £;jAx"^/^cos$ Q 
V2fro {rl + J^f 99 



{S^ + 1 - cos 9) 



-1/2 q^ErVoX ^ — i 



2 {rl + J2) 



Po 



(84) 



-3/2 



the limit A — and the first term on the 



As we have seen before, by Eq. ^ {5"^ + 1- cos 9) 
right hand side will vanish. The second term will also give no contribution to the regularization parameters because 
(x "^^^ cos <i>^ = 0. Only the last term, which is ~ will give non-zero contribution according to the argument 

in the analysis presented above (see Eq. ([7^ ). Using Eq. (HTl) in the limit A ^ and taking "()" process, Eq. 
becomes 



o{l + JVrlf^U^ 



(85) 



The identity {x 



2F1 {p, i; 1, a) = Fp, with a = P/ (r^ + J^) is taken from Appendix C 



of Paper I (author?) [7], and we take the limit 9^0 



^un Q^(i)[e 1] 



1 g2 ErF y2 
2r 



o (1 + jyri) 



3/2 ■ 



Now the remaining part is 



3q2 



Po' [dt{p')]Vuil^, 



(86) 



(87) 



where 



[dt{p')]Viii\ 



2ErA 
1 



Mb? 



f 2.P 



5M\ JrA^ 



(0- 



roAl 



-roJf {(f> - (l)of ~ ToJr {(j> - (j)o) (o - 

Taking similar procedures as above, the non-vanishing contributions turn out to be 

3 



+ O[{x-xo)% 



\a->o 



Hm -q'^EJ^rroPr, ^ cos^ <i> sin"' 9 
lim^^-^!S.(x- 



A^o 2 To 1 + J^/rl 
3 g2 Er I 



2r, 



o (1 + JVP) 



2^,3/2 



1 + J^lrl 



where all other terms than ^ p^^ again have been dropped off during the procedure since they vanish either 
limit A — >■ or through the "()" process. Then, using the identity = 2F1 (p, ^; 1, a) = Fp, we have 



3g2 



Er 



By combining Eqs. (|55)) and (PU)) . we finally obtain 



5/2 



Bt = ^Er 



3/2 



3F. 



5/2 



(l + JVr2)3/2 2(l + J2/r2) 



5/2 



3. Br -term: 



From Eq. we start with 



~Pl 



^ 2 pI 4 

Then, following the same steps as taken for the case of Sfterm above, we obtain 



Br — 

r 



1/2 



(1 - 2/-lf2)i^3/2 , 3/-lf2F5/2 



(l + J2/r2)^/2 2(l + J7r2)3/2 2(l + j2/r2) 



2^5/2 



5. Bjy-term: 



Again, from Eq. (|64p 



Then, similarly we can derive 



Br/: 



t = to 



T 



Fl/2 — F3/2 3(^5/2 — F3/2) 



(l + J2/r2)i/2 2(1 + J2/r2)3/2 
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4- Bg-term: 

As Ag vanishes, so should Bg. From 

]-g|-2 p3 +4 1' (96) 

one finds that there is no term like ^ p~^^: all terms are either like ^ /\2n^~2ri+i j^jj^g ^ A^"~^ sin cos $7/5^"+^ 
(ri = 1, 2), which vanish in the limit A — > or through the "()" process. Thus 

Be = 0. (97) 



C. Ca-terms 

We have mentioned before that Ca-terms, which originate from e°-term in Eq. (|43|1 . always vanish. This can be 
proved by analyzing the structure of e^'-term. First we specify the e^-order term ^ot daiX / p)\^^^ ™^ ^ Laurent series 
expansion and define 

^"t^ ^ " ' 1~2^^ + 4 16 Jl 1 ■ 

Generically, this can be written as 

3 2n+l[fc/2] A2ri+l-fc/i ^ ( 

Qai^"] - E E E — (99) 

n=l k=0 p=0 P° 

where A = r — Tq, and Cnkp{a) is the coefficient of each individual term that depends on n, k and p as well as a, with 
a dimension TZ''"^ for a = t, r and T?,*^"^ for a = 6, (p. 

The behavior of Qa[e°], according to the powers of each factor on the right hand side of Eq. ((99)) . is 



Qa[e"] - p-(2"+i)A2"+i-^ (0 - <^o)'=-2^ - f ) ^^ (100) 



where s = k — 2 for a = t, r and s = fc — 1 for a — 9, (f). Following the same procedure as in the beginning of 
Subsection IV B[ Eq. (llOOp becomes 

ga[e"] ^ p-(2«+i)^2n+i-2p-. (sine)'^'+' (sin$)'^ (cos$)^7^^ (101) 

where a binomial expansion over the index i = 0, 1, • • • , fc — 2p is assumed, and s — 2p + i — 2 for a = t, r and 
s = 2p + j — 1 for a — 9, (p. Here we have disregarded any by-products like 0[{x — Xo)'"'^^^^^] and 0[{x — Xo)'^^^'^], 
which originate from (0 — 0o)'^ and (6' — ^) ^, respectively when we rotate the angles: by putting them back into 
Eq. (jlOOp we simply obtain e^-terms, which would correspond to 0(£~^) in Eq. PU)) and should vanish when summed 
over £ in our final self- force calculation by Eq. © . Then, the rest of the argument is developed in the same way as 
in the beginning of Subsection IV Bl 

(i) * = 2j + 1 (j = 0, 1, 2, • • •) 

The integrand for "()" process, F($) = (cos$)^^^^ (sin<I>)^'' has the property F($ + tt) = -F{^)- Thus 

{Qa[e°])^0, (102) 



(ii) ^ = 2j [j = 0, 1, 2, • • •) 
We have 



Qa[e"] (sin$)^P(cOS$)^^'p-2(n-g)-1^2(„-g) + l-^s^ (^Qg) 



where q = 0, 1, • • • , p + j is the index for a binomial expansion and s — —2 for a = t, r and s = —1 for a ~ 9, (j). 
Here we can guarantee that n — q > — i, i.e. n — q = 0, 1, 2, • • • since < q < p + j — p + ■^i, < i < k — 2p, 
and p < k < 2n + 1. Then, Eq. (|103p can be subcategorized into the following two cases; 
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(ii-1) n — q > I 

By Eqs. gl|), and (USD 



]a[e°] - (sin$)'^(cos$)'^ A2p^(cose)7^" ^0, 



(104) 



(ii-2) n~ q^O 

By Eqs. dSSD, gl|), and dEj 



Qa[e°] - (sin$)^P {cos^f^ AP^(cose)7^" ^ 0, 



(105) 



where s — —2 for a ^ t, r and s = — 1 for a = 9, (/). 

Clearly, in any cases the quantity Qaie*'] does not survive, therefore we can conclude that Ca-terms are always zero. Q. 
E. D. 

Also, this justifies the argument that we need not clarify the term 0[{x — Xo)^] in Eq. (|2ip and its contribution 
to which is 0[{x - Xo)^] in Eqs. ^ and ^ in Section |IV] or Vw in Eqs. ^ and dH]) in Section El by the 
analysis of the generic structure given above, — i daViyl^^^^ /pi or | [da (p^)] 'Piv\^_^ / Po would simply vanish in 
the coincidence limit gardless of what Viy is. 



VI. DISCUSSION 



The scalar field in the flat spacetime limit 

One interesting fact is that when we consider flat spacetime as the background our singular source field -0® may 
be determined in a completely different way from the case of curved spacetime. This flat spacetime version of V'^, 
though obtained using a different method, should agree with its curved spacetime version when the flat spacetime 
limit is taken. 

Without introducing the special coordinate frame like that of the THZ coordinates as employed in Section IIVI this 
argument can be shown rather straightforward: the retarded field, which is equivalent to the singular source field in 
flat spacetime, can be computed directly using the coulomb potential and considering the special relativity. Here we 
will focus on p — \x ~ Xo\ between the field point x and the source point Xo- This will be expressed in terms of the 
Cartesian coordinates first. Then, we will switch from the Cartesian coordinates to the spherical polar coordinates 
and arrange the terms in a polynomial according to their orders, where quartic or higher order terms will be separated 
as errors. Finally, the expression of p^ for ip^ct in flat spacetime obtained thus will be shown to agree with the flat 
spacetime limit of p^ for t/^g in Eq. (|29p . 

In the frame of reference in which the particle is always at rest at the point x^, the retarded field at x is 

, Q 

Aet = 7^ 



We boost this frame such that in the new frame of reference the particle is moving with the 3-dim velocity /3 (see 
Ref. (author?) [HI): 



jt = X ~ Xo 



where = x' — xi and T' = t' — t'„. Then, we have 

-7-1 



p'=\x-Xo\={{x'-x',)+P 



I3.{x'-x'j ~l[t'^Q 
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Dropping the ' notation and expanding the terms inside { } out 



Now we need convert this expression into the spherical polar representation. First, using the cosine law, the distance 
between x — (r, 9, 0) and Xq — (ro, ^o, (t^o) can be expressed as 



where 



If — Xo|^ = +rl ~ 2rro cost}, 
cos ■& = sin 6 sin 9^ cos (0 — 0o) + cos cos 9^, 



which is obvious from the trigonometric rule. In particular, for the particle moving along an equatorial orbit {9o = 7r/2) 
we may rewrite the above as 

\x-Xof = {r - ro)^ + To (r - To) (j) - + ro {r - ro) {(f) ~ (j)of + (9 - 

+rl {<j>-cj,of + 0[{9-TT/2,cl,-cj,o)% 

where all the trigonometric functions of the small arguments are expanded in Taylor series up to the cubic order. 
To compute f3 ■ {x — Xo), first of all one need change the basis vectors from {x, y, z} to {fo, 9o, 4'o} via 

cos 0o — sin 0o \ / 
sin0o cos(/)o 9o 

Then, be rewritten as 

X — Xo — (r sin 9 cos 0) o cos (j)o — 0o sin (f>o^ + {r sin 9 sin cf)) □ sin 0o + 0o cos 



-r cost 



' o ■ 



Also, in the new basis the 3-dim velocity of the particle moving in the equatorial plane is expressed as 

P = Pr^o + P<P<i>o- 

Thus, taking a dot product of /3 and x — Xo gives 

(3 ■ (x — Xo) — I3r[r sin 9 cos {4> — (j)o) — ro] + /?0r sin 9 sin {4> — (j)o) ■ 

Again, taking a series expansion of this quantity around the source point Xo — {ro, f , (/>o), sufficiently up to the 
quadratic order, we have 



(3 ■ (x - Xo) = Pr {r - To) + /30ro ((/> - (t>o) + {r - ro) {(t> - 0, 

-\rolir (0 - 0o)' + 0[{9 - - 4>o)% 
Finally, putting all the terms together we obtain the following 



p 



{E^ -l){t- to)' - 2Er [t - to) {r ~ r,) - 2EJ [t - Q (0 - 0o) 

+ (1 + r') {r - ro)' + 2Jr (r - ro) (0 - 0o) + {rl + j') (0 - cpof + - f )' 

(t - to) {r - ro) {(t> - M + roEr (t - to) (0 - 0o)' + ToEr (t - to) - 



2EJ 



2Jr 



2J2 



+ (r- - ro)' (0 - 0o) + ro 1 - r' + ^ (r - ro) (0 - 0o) 



+ro(l-f')(r-ro) [9-^)' 



— roJr ( 



(/.o)'-roJr(^-0o)(0-|)' 



+0[(t - to, r - ro, - tt/2, - (/)o)4]. 

This is exactly equal to the flat spacetime limit of for tp^ when M = in Eq. ([29l) . with appropriate replacement 
for the coefficients, using E = (dt/dr)^ = 7, J = j-q {d(t)/dT)^ — jro(3^, and r = (dr/dr)^ — jfir together with the 
identity P = E^ - f {l + .P/rl). 
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Appendix A: THE THZ NORMAL COORDINATES 

In Section ITVl we have introduced the Thorne-Hartle-Zhang's normal coordinates to simplify the description of 
as if it was measured by an observer who travels on a particle moving in curved spacetime: in this coordinate system 
{X^}{A = 0, 1, 2, 3), = q/p + 0{p^/TZ^), where = SjjX'X^ {I, J = 1, 2, 3). As presented by Eq. the 
specification of 0[{x — Xq)^] in p^ (thus of 0[{x — Xo)^] in X^ in Eq. (PTjl ) is not necessary for our current mode-sum 
regularization scheme: it proves not to contribute to the regularization parameters Aa, Ba, and Ca- 
in the practical calculation of self- force, however, it will be very useful to extend our scheme to the next orders, 
i.e. to Da or higher terms. If we extend Eq. ()43|) . say, to e^-term, it would generate the next-order regularization 
terms -2^Da/[(2£ - 1)(2£ -I- 3)] in the place of 0{e-'^) in Eq. Strictly, these terms would give non-vanishing 

contributions to the self-force since the sum oi -2V2Da/[{2£ - l){2e + 3)] is taken over many but finite number of ^'s 
in actual numerical calculations. 

In order to determine Da or higher terms, the knowledge of 0[{x — Xq)^] in p^ and thus of 0[{x — Xo)*] in X^ will 
be required. As it should be an important tool for this purpose, we present below a more detailed description of the 
THZ normal coordinates = 0, 1, 2, 3) for a particle moving along a general orbit (confined to the equatorial 

plane for convenience) about a Schwarzschild black hole [22j: the expressions are given in terms of the Schwarzschild 
coordinates x°' — (t, r, 9, (f) and specified up to the quortic order 



T=X^ = -u,A [X^ + a^BcnX^X^X^ + bcdeX^ X^] + 0{X^), 
X' = n^P^ [X^ + k^bcdX^X^X^ + \^ bcdeX^ X^ X^ X^] + 0{X^) 



(Al) 



with 



X^ = M^aix'' " + \M^a " - O + 0[{x - x,)% (A2) 

a^BCD ^ ^r^Q,«|j7r^B7r«c^«z, + 3^'3B^VD 

+3^«b/i^c/i'^d + h^'Bh^^ch'^D) , (A3) 

GTr'^BTr'^ch'^Dh^E + Aw"^ Bh"" ch"^ nh"^ e + h"" Bh'^ ch"^ oh^ e) 
5 

"168 



R^PQR,s\^TTQ^hPBh''chDE. (A4) 



D 



n^BCD ^ ir^Q,^|j7r^B7r«c^^D + 3^QB^W 

+ iir''Bh^ch'^D + hPBh^ch^D) 

- R^'pQbL { b^ETT^'^h^BhcD + Is^BTr'^Eh^ch'^D] , (A5) 



PQfllo ET^ II Bll-CD^-^^ 

X^BCDE = ^ r^Q^^5U7r^B7rQc7r^D7r^ij + 47r'3s7r^c7r'^D/i^£; 

Qn'^B^'^ch'^Dh'^E + An'^Bh^'ch'^Dh'^E + h'' Bh'^ ch'^oh^ e) 

+ ±S^pn^Qh^Bh'chDE + ^S^BTr'^Fh^ch^Dh'^E + pTT^' bH^ choE 

24 24 o3 



(A6) 
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where A, ■ ■ ■ , E a.nd P, ■ ■ ■ , S = 0, 1, 2, 3, and /, J, K, L = 1, 2, 3. For Eq. dJT]) we have 

< - (/-^/'i?, f-'^'r, o) , (A7) 



n. 



n(2)^ = (_:^ 1 + — ol fA9) 

1^ ro'ro(i? + /V2)''+,2(;-i/2^+l)'^j' 

n(3)^ = (0, 0, 0, 1) , (AlO) 

where f = (^l ~ V^)' ^'^^ ^ = = (1 ~ 2A//ro) (dt/dr)^ (r: proper time) and J = = {d(j)/dT)^ are the 
conserved energy and angular momentum in the background, respectively, and r = — {dr/dr)^. In Eq. (|A2[) we 
define 



Co = (io, ro, |, 0o) , (AH) 



M^, = diag _^„J , (A12) 

along with the non-zero Christoffel symbols at Xo in the Schwarzschild background 



^t'-lo ~ 712' -"""lo - "2"' -'"^'-lo - 7~2' ^ee\o - /^o' ^4"}>\o ^ ^'^re\o - ^t4> - ~- (A13) 

J ' o ' o J ' o o ; n 



The quantities ^|^, R^bcd\^, ^bc de\o R'^bcd,e\^ in Eqs. (|A3P - (|A6p are evaluated from the initial static 
normal coordinates {^"*} represented by Eq. (jA2p . They follow the identities 

TbCdIo ^ H^BCD + H^cBD — Hbc^d, (A14) 

R^BCd\^ = Hbc^ D — H^CBD — Hbd'^C + BiBC^ (A15) 

TbCDbIo ~ 'i {H^ BCDE + H^CBDE — Hbc^ De) , (Al6) 

^'^BCD,£;|q = "i {Hbc^ DE — H^CBDE — Hbd'^ce + bibce) ■, (A17) 

where the building blocks Habcd and Habcde are taken from — r/^is + HabcdX'-^ + HabcdeX'^' 

(the linearized gravity in the geometry of {X"*}) and have symmetric properties Habcd ~ H(^ab){cd) and Habcde = 

H(ab)(cde)- The non-zero Habcd and Habcde turn out to be 



-"0000 — ^TT' ^0011 — T — q 3~ ' -"0101 — "TT' ^0202 — -"0303 



1 (2M 3M2 



_ Af2 /2M A'P\ _ _ f _ _ I ( I M 

-ffllOO — TT-> -^1111 — T — 5 T ■ -^1122 — -ffll33 — -^1212 — -H^1313 — ' 



and 



2 _ 1 /2M 

_ / 1 _ / 

-^2222 = ^, -^2233 = -^2323 ^ TTT ) ^^3333 



1 /2Af2 3Af3\ 1 /2M SAf^ 

I — I r7 nm 1 1 I - — -I— 



-ffoOOOl = — „ „ /o — P 3— , -ffool 



3J3/2 ^ ^5 ^6 ; ' _^3/2 ^ ^4 ^5 ^6 ^ ' 

1 /2A^ 3A4"2\ _ Ai"3 _ 1 /4Af2 3Af3\ 

-ff00122 - -ff00133 - (^^;4 ^ j , -^01000 - " j57i;|' ™ll - -gy^ (^^:5 j ' 
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01022 



A/2 



01033 — 



3/1/2^5 



02012 



03013 — 



6/1/2 



M 3AP 

^4 



Hi 



1 /2A/2 3Af-^ 



1001 



,i?i 



1 /2A/ GM^ 



^^11122 — -fflll33 — 



3/3/2 y ^5 ^6 y ^^^^^ J3/2 y ^4 

1 /4 14Af 15Af2' 



3/1/2 ^^3 ^4 



i?l 



1 /Af Af2 



2002 



6/1/2 ^sy 



^^12 



6/1/2 



4 IIM 9M2\ 

::3 - -IT- + -;:5- ) ' ^12222 



/ 



1/2 



A I : -"12233 



pl/2 



Af 



H 



1 /Af Af 



13003 



//-2 



6/1/2 y^4 ^5 



1 / 4 llAf 9Af2 



13113 — 



6/1/2 y^3 ^4 



pl/2 



H 



13223 



J ] , -"13333 — 



/1/2 / 1 



M 
74 



22122 



2/1/2 / 3^,j 



33133 



2/ 



1/2 



22133 



3r3 



-^23123 — 



/1/2 / 1 



3A/ 



According to Ref. ((author?) [13), one can show the foUowmg out of the two geometries, the THZ {A""^} and the 
Schwarzschild {x"} 



9a;" dx'' 



--45 .-0(^-2), 



(A18) 



where tilde denotes the THZ geometry and 0[{x — Xo)^] is converted to 0{X^) via the inverse transformation of 
Eq. (jA2[) . The metric perturbations 0{X^) in the last line of Eq. (jAlSP should take specific forms to satisfy the 
properties of the THZ coordinates as mentioned in Section IIVI Our results show that 



500 



— 1 — £klX X 



K 



^ c vK vL -wM 



-^e-iKpB^ lX^ — —EklX^X^Xi + —p^EkiX 



K 



1 



eiKpE'^LuX'^X'^X^' + 0(p77^"), 



K vL 



gij = Sij — 5ij£klX X 



f-iKpS^ lX^ X^Xj 
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-l^p^KPiBj^^X'' - UijEkluX'^X'^X^' + 0{p^lTl^), 



(A19) 



(A20) 



(A21) 



where = X^ + y"^ + Z"^ and indices /, J, K, L, M, P = 1, 2, 3. The external multipole moments are spatial, 
symmetric, tracefree tensors and are related to the Riemann tensor evaluated on the particle's worldline by 



1 



Bij 

ElJK 
BjJK 



2^1^^ RpQJO 



Vif i?o/0.7 



STF 



3 r ~ nSTF 

- tl^'^V K RpQJO 



(A22) 
(A23) 
(A24) 
(A25) 
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where STF means to take the symmetric and tracefree part with respect to the spatial indices /, J, K. The dot 
denotes differentiation of the multipole moment with respect to T along the particle's worldline. One should see that 
£ij ~ Bij ^ 0(1/7?.^) and £ijk ^ Bjjk ^ £ij ^ Bu ^ 0{1/'R?) for consistency of the dimensions. The fact that 
all of the above external multipole moments are tracefree follows from the assumption that the background geometry 
is a vacuum solution of the Einstein equations. These results agree with Eqs. (17) and (18) of Paper I (author?) |3| 
or Eqs. (3.26a)-(3.26c) of Zhang (author?) [9] to the lowest a few orders. 

Appendix B: HYPERGEOMETRIC FUNCTIONS AND REPRESENTATIONS OF REGULARIZATION 

PARAMETERS 



In Section El we define 



X = 1 - asin^ $ (Bl) 



with 



J2 



r2 + J2 • 



(B2) 



And we use 



TT 



-it/2 



fl-asin^$) 







2 

In particular, for the cases p = ^ and p ^ — \ we have the following representations 



2FAp,]--l,a]=Fp. (B3) 



and 



J^i/2=2fi =5i^(a) (B4) 



= 2^1 ( i l;a) = ;^^(«), (B5) 



where K{a) and E{a) are called complete elliptic integrals of the first and second kinds, respectively. 
If we take the derivative of F1/2 with respect to fc = ^Ja via Eq. (jBSp . we obtain 

dk k k ' 



or using Eq. (jB4 



However, Ref. (author?) [1^ shows that 



dK K ttF. 



3/2 



dk k 2 k 



dK _ E K 
'dk ~ fc(l - fc2) ~ 

Thus, by comparing Eq. (|B7[) and Eq. (jBSp we find the representation 

2 E 2 E 



(B6) 



(B7) 



(B8) 



^3/2 = = -1 • (B9) 
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Further, we can also find the representation for First, taking the derivative of ^3/2 with respect to fc = ^/a 

via Eq. (jBSP gives 



3/2 



dk 



(BIO) 



Also, using Eq. (|B9p together with Eqs. (jB3P - (jB5p . another expression for the same derivative is obtained solely in 
terms of complete elliptic integrals 



3i^3/2 _ 2 (1 + ^ - (1 - k"^) k 



dk TT 

Then, by Eqs. (HH, (|BT0)) . and, (|BTT|) we find 

-P5/2 = ^ 



fc(l - 



2{2-a)E K 



{1-ay 1 



(Bll) 



(B12) 



Now, using Eqs. (jB4[) . (|B9p . and (jB12p . we may rewrite the non-zero 5-terms of regularization parameters, Eqs. 
CT - p^ in Section Iml as 



-,2 Ef 



K{a) - 2E{a) 



(B13) 



Br = 



{P - 2E^) k{a) + {P + S^) E{a) 
^ TT{l-2M/r,){l + J^/rl) 



(B14) 



B, 



2 r K{a) - E(a) 



^ ro[J/ro){l + JVrlf'^' 
which are exactly the same to the results of Barack and Ori (author?) [lol |. 
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Or alternatively, one can use the arguement ^ |^ P^(cos O) = as B — >■ 0, to show that this part does not survive at the 
end. ^ 

The by-product from 0[{x — Xo)'^] in Eq. (|53p . that is to say, — ^ p~'^0[{x ~ Xo)'^] is suppressed since this can be categorized 
into the Ca-term group (see Subsection IV C|) . 
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[20] The by-product from 0[{x — Xo)^] in Eq. (|57ll . that is to say, — ^ p~'^ 0[{x — Xo)^] is suppressed since this can be categorized 
into the Ca-term group (see Subsection IV C[) . 

[21] Developing our mode-sum regularization scheme further, we may include e^-term in Eq. (|43p and this would generate the 
next-order regularization terms —2y/2Da/[{2£— 1)(2£-|-3)[ in the place of 0{£~'^) in Eq. (|10p . Rigorously, their contributions 
to the self-force would be non-vanishing since we are taking the sum of -2V2Da/[{2£ - 1){2£ + 3)| over many but finite 
number of i's. Thorough discussions on D^-terms for a circular and for a general orbit cases are found in Ref. (author?) 
Q and Ref. (author?) 17![, respectively. 

[22[ For derivations of the THZ-coordinates and for discussions on their application to Da-terms, readers may be referred to 
Ref. (author?) 7[ and Ref. (author?) [17[ for a circular and for a general orbit cases, respectively. 



